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Abstract 



The Crab Nebula is one of the brightest and most stable sources in the X-ray sky. Year-scale flux 
variation from the object was recently revealed in the hard X-ray band by four satellites. This marked 
the first detection of year-scale variability from pulsar wind nebulae in the hard X-ray band. The Crab 
Nebula has been observed at least once a year for calibration purposes with the Suzaku Hard X-ray 
Detector (HXD) since its launch in 2005. In order to investigate possible spectral changes as well as flux 
variation, the archival data of the HXD were analyzed. The flux variation reported by other instruments 
was confirmed in the 25 - 100 keV band by the HXD in a few percent level, but flux above 100 keV did 
not follow the trend in variation below 100 keV. The hardness ratios produced utilizing the PIN and GSO 
sensors installed in the HXD exhibit significant scattering, thereby indicating spectral variations in the 
hard X-ray. The spectral changes arc quantified by spectral fitting with a broken power-law model. The 
difference between the two photon indexes of the broken power-law model in harder and softer energy bands 
f*^) ■ is in the range of < 2.54. Taking into account flux variation of 6.3% and spectral variation time-scale of 

a few days, multi components of the broken power-law-shaped synchrotron emission with different cooling 
t ■ times are suggested. 

' Key words: radiation mechanisms: non-thermal — ISM: individual (Crab Nebula) — pulsars: indi- 

O ! vidual (PSR B0531+21) — X-rays: stars 

CO 

1. Introduction ages. The flux of the pulsar component accounts for ^20% 

of the total flux with most of the rest stemming from the 
' The Crab Nebula is one of the most famous pulsar wind PWN in the hard X-ray band. 
*_h . nebulae (PWNe) which is located in the center of a histor- The pulsar essentially drives the electromagnetic radia- 
ical supernova first recorded in Japan and China in 1054 tion as follows: energetic flows of electrons and positrons 
(Duyvendak 1942). An energetic, bright pulsar (known from the pulsar (pulsar winds) reach and interact with 
as Crab pulsar=PSR B0531+21) exists in the center of interstellar mediums to induce a termination shock. The 
the Crab Nebula. The Crab Nebula has been observed in high energy particles are believed to radiate synchrotron 
the entire observable electromagnetic wavelength from ra- emission outside of the shock. This radiation is called as 
dio to TeV gamma-rays. From the non-thermal spectrum PWNe. Since the origins of the emissions of pulsars and 
and strong polarization (Wcisskopf ct al. 1976; Forot et al. PWNe are the rotation of magnetosphere of the neutron 
2008), emissions from the object are interpreted as syn- star, they had been believed to be stable in principle, 
chrotron radiation in the X-ray band. Emissions from the Data stored in ~12 years recently indicated 10 percent 
direction of the object are divided into PWN and pulsar flux variation in 10 - 300 keV with hard X-ray instruments 
components. These two components can be separated not (Wilson-Hodge et al. 2011). As reported by Wilson-Hodge 
only by imaging analyses but also by phase-resolved tim- et al. (2011), the pulsed flux of the Crab pulsar in the 2 
ing analyses if we assume the emission from the pulsar is - 100 keV band slightly decreased by 0.2% yr _1 , which is 
almost the same as the pulsed emission. Therefore, they negligible for hard X-ray variation. Therefore the varia- 
are separately observed and discussed by many authors tion seen in the Crab observations is supposed to be due 
even in the hard X-ray band, where it is hard to get im- to that of the flux of PWN. Currently, there are no rea- 
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sonable models to account for the variation of the flux of 
PWN. Wilson-Hodge et al. (2011) also suggest the soft- 
ening phenomena in the light curves of < 100 kcV with 
Fermi/GBM and INTEGRAL /ISGRI. 

As indicated by Jourdain & Roques (2009) and Ling 
& Wheaton (2003), the X-ray spectrum from the Crab 
Nebula is described by the broken power-law model, whose 
break energy is at about 100 keV. This break can be 
explained by electron cooling via synchrotron radiation. 
Thus the break point and the index may be changed in 
accordance to the flux change. In fact, Ling & Wheaton 
(2003) reported variable broken power-law spectra with 
CGRO/BATSE, while Jourdain & Roques (2009) could 
not confirm the phenomenon. Thus, the existence of the 
variation in the broken power law spectrum was unclear. 
Here we analyzed the archival data of the Crab Nebula ob- 
tained by the Suzaku Hard X-ray Detector (HXD), which 
offers the best sensitivity in the range of 10 - 600 keV. 

Since the promising variability could be so small com- 
pared with the systematic errors in the standard calibra- 
tion of the instruments, in principle, we have to clarify 
the validity of the energy response matrix including the 
effective area before the spectral analyses. This paper 
reports on spectral variations with close attention paid 
to instrumental calibration, and is organized as follows: 
Section 2 discusses the observations and the method of 
data reduction employed. Section 3 describes the results 
obtained. In Section 3.2, we first examined the spectral 
changes without model fittings in order to reduce possi- 
ble errors from the uncertainty of the response matrices 
and the method of model fittings. Then in Sections 3.3 to 
3.5, we performed model fittings to the spectra in order 
to quantitatc the spectral changes. Section 4 describes 
the analysis and results regarding the pulsar component. 
Finally, Section 5 discusses the origin of the observed vari- 
ability. 

2. Observation and Data Reduction 

2.1. Observation 

We used data observed by the PIN (10 - 70 keV) 
and GSO (Gd 2 Si0 5 (Ce); 40 - 600 kcV) of the HXD 
(Takahashi et al. 2007; Kokubun et al. 2007) on board 
Suzaku (Mitsuda et al. 2007). The HXD has two ad- 
vantages in investigating hard X-ray spectral variability: 
wide-band sensitivity in 10 - 600 keV, and small back- 
ground rate and high reproducibility (i.e., smaller system- 
atic errors) of the background model. Suzaku also car- 
ries CCDs called the X-ray Imaging Spectrometer (XIS; 
Koyama et al. 2007). However these devices are difficult 
to use for this study due to a pileup problem when ob- 
serving the bright object like the Crab Nebula. For this 
reason, we did not use the XIS data in this work. 

There are two typical observation attitudes for Suzaku 
corresponding to the optical axes of the XIS and HXD 
(respectively called "XIS nominal" and "HXD nominal" ) . 
Suzaku observed the Crab Nebula with both nominal po- 
sitions to check the effective areas and energy responses. 
Observations for this purpose were made 17 times from 



July 2005 to March 2012, and all data were public soon 
after the observations. In order to avoid relative large un- 
certainty in calibration, we excluded the datasets other 
than nominal pointing positions; i.e., offset observations 
of the Crab Nebula, which is used to define the optical 
axes at the initial phase of the satellite. 

The effective area of the HXD is different by the nom- 
inal positions simply due to the transparency of its fine 
collimator. Takahashi et al. (2007) and Kokubun et al. 
(2007) show that the ratio between effective areas of the 
HXD and XIS nominal positions is about 92%, and the 
differences in energy dependencies of the effective area be- 
tween them is within 1% in 15 - 70 keV and 2% in 50 - 
600 keV. Every response matrix is defined for a period, 
"epoch", which is divided by the change of operational 
parameters; combination of a bias voltage and a set of 
lower discriminator (LD) level. Table 3 of Nishino et al. 
(2010) summarizes the history of these operations. 

2.2. Data Reduction 

The gain changes of the PIN are within 1% (Nishino et 
al. 2010). We also confirmed the PIN gain stability with 
accuracy of < 0.7% until 2012 September after the study 
by Nishino et al (2010), by measuring the Gd K lines 
produced in GSO crystals which are irradiated by X-rays 
from objects (Kouzu 2013). The gain of the photomul- 
tipliers for GSO changes in various timescales, so that 
the GSO gain history files are constantly revised along 
with calibration results. The energy scale in < 100 keV 
was also improved by the FTOOLS hxdpi (Yamada et al. 
2011). All HXD data were reprocessed with aepipeline 
version 1.0.1 in the HE AD AS 6.11 software package using 
CALDB 2011-09-15; therefore, the GSO data were applied 
to hxdpi version 2010-01-12 and the gain history gsogpt 
in this work. 

Background for the HXD data is generated as a 
synthetic model that accounts for the time-variable 
particle background (the "non X-ray background", 
NXB; Fukazawa et al. 2009). Specifically, we 
used the "tuned" PIN NXB models of version 
2.0 (METH0D=LCFITDT , METH0DV=2.0pin0804). We 
also used the GSO NXB models of version 2.5 
(METH0D=LCFIT , METH0DV=2.5ver0912-64) for observa- 
tions before November 2011, and those of ver- 
sion 2.6 (METH0D=LCFIT , METH0DV=2.6ver 11 10-64) after 
November 2011. The uncertainty of the NXB models of 
the PIN is within 3% (Fukazawa et al. 2009), which is neg- 
ligible as the Crab source count rate is 10 times greater 
than that of NXB, even if at 70 keV. The uncertainty of 
the GSO NXB is less than 1% (Fukazawa et al. 2009). 

In the estimation of the flux from the object, it is impor- 
tant to calculate the dead time (or live time) of the obser- 
vation. The dead time of the HXD-PIN and GSO is cal- 
culated regarding the following three processes: (1) event 
processing stopped by another process running previously 
triggered, (2) an event data packet discarded in commu- 
nication between the HXD analog electronics and digital 
electronics, and (3) space packets containing events dis- 
carded due to limited bandpass between the digital elec- 
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tronics and data processor. In particular, the HXD team 
tunes the parameters of analog electronics, and so the case 

(2) only occurs just after South Atlantic Anomaly pas- 
sages. Also, we discard events during the periods in case 

(3) by the FTOOL, hxdgtigen. The dead time is calcu- 
lated and stored with the pulse height data obtained by 
the FTOOL hxddtcor. The dead time fraction is deter- 
mined with accuracy of < 0.2% at an observation duration 
of > 10 ks (Kokubun et al. 2007). 

3. Analysis of Total Emission 

3.1. Flux Variation 

In order to examine the time variation of total flux 
(pulsar + nebula) reducing effects by any systematic un- 
certainty of the instrumental responses, we evaluate the 
count rates of individual observations by subtracting both 
the non X-ray background and cosmic X-ray background 
from events. We also have to note that different bias 
volatages for PIN were applied during epoch 1, epoch 2 
and after epoch 3, as represented in Table 3 of Nishino 
et al. (2010). Effective areas for each epoch is calculated 
based on ground tests (Sugiho et al. 2001; Nishino et al. 
2010). In order to compare the measured flux in each 
epoch, here we normalized the photon flux using the cal- 
culated effective area for each epoch. In addition to the 
effect of bias voltage, the effective area decreases at < 
25 keV because of LD of each PIN detector. Besides, the 
detector's responses above 55 keV were affected by the 
degraded depth of depletion layers, which is adjusted by 
using the calibration data-sets on the Crab Nebula itself. 
Therefore, we only used the 25 - 55 keV band where the 
effective area of the response matrices changes within 1% 
for observations after PIN epoch 3. Responses of this en- 
ergy band depend on neither signal processing nor time 
degradation, and are only defined simply by the physical 
process; cross section of the photo absorption process, or a 
stopping power of Si of the PIN diode (Nishino et al. 2010; 
Kokubun et al. 2007). The count rates are normalized by 
the effective area calculated for each observation, taking 
into account the difference in effective area according to 
different nominal positions. 

Figure 1 shows the variation of total count rates with 
the HXD and other X-ray detectors normalized with the 
count rates in MJD 55000 - 55500. The count rates in 25 
- 55 keV and that in 50 - 100 keV are consistent with the 
year-scale trends by other satellites except for the data 
after MJD 55500. On the other hand, the trend of the 
GSO data in 100 - 500 keV deviates from trends of the 
other energy bands even when the data include statistical 
errors (la) and systematic errors (1% of the NXB). In 
addition to the year-scale changes, a day-scale variation of 
~1% with the PIN and GSO is also represented, which is 
derived from a ratio between count rates of 2006.03.30 and 
2006.04.04 whose configurations (PIN epochs and nominal 
positions) are the same (see Table 1 for the time intervals 
and count rates). 
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3.2. Hardness Ratios 

To clarify the year-scale and day-scale variations 
(Section 3.1) without model fittings, we produced the 
hardness ratios of the PIN and GSO datascts. Figure 2 
shows the count rates vs. hardness ratios of the PIN and 
GSO. The count rates (horizontal axes) are normalized 
by the effective area calculated for each observation, tak- 
ing into account for different nominal positions. To avoid 
mixing the data of different detectors (i.e, different energy 
response matrices), here we plotted hardness ratios of the 
PIN and GSO individually. When observing such bright 
sources as the Crab Nebula, Gd K lines (~43 keV) from 
the GSO had not been reproduced with sufficient accu- 
racy in the response (Nishino et al. 2010). We did not use 
the energy range of 40 - 45 keV to avoid the Gd K line 
structure. 

In order to check correlations between hardness ratios 
and count rates, we fitted the data with linear functions, 
and then got x 2 /d.o.f = 16.72/10 (PIN) and 47.96/15 
(GSO). Hypotheses of linear correlations are rejected on 
91.9% (PIN) and > 99.9% (GSO) confidence levels. In 
other words, both the hardness ratios and individual count 
rates vary significantly, but neither shows any unique cor- 
relation. Moreover, the hardness ratios of the soft band 
and hard band behave differently, suggesting a variable 
break energy or photon indexes in spectra which they are 
described by the broken power law model, as indicated by 
Ling & Wheaton (2003). 

3.3. Individual Spectra 

The level of the observed variation (~ 6.3%) of PIN 25 
-55 keV count rates, as evaluated in section 3.1, is fairly 
close to the sum of official values of the systematic er- 
rors of released response matrices (gain stability), NXB, 
and dead time correction (< 1.2% total at < 100 keV; 
Section 2.2). Here we summarize evidences showing that 
the above analysis is reliable. As the response matrices of 
the HXD-PIN (25 - 55 keV) and GSO are constructed only 
in line with Geant4 simulations with parameters obtained 
in ground tests (Tcrada et al. 2005), they have not been 
tuned or calibrated with the Crab Nebula. Thus we rea- 
sonably consider the response matrices to be independent 
of the flux variation of the Crab Nebula. As mentioned in 
Section 3.1, the HXD flux trend was consistent with other 
satellites at 3±1% (PIN 25 - 55 keV), 1±1% (GSO 50 - 
100 keV), and 4±1% (GSO 100 - 500 keV), taking into 
account day-scale variability as the second errors (Fig. 1). 
These facts support the conclusions that effective areas 
of the PIN 25 - 55 keV, GSO 50 - 100 keV and 100 - 
500 keV are essentially stable within 4%, 2% and 5% at 
most respectively, and that the relative flux variation is 
real. 

Before evaluating individual spectra, we first performed 
a model fitting of the averaged spectrum of all HXD ob- 
servations in 25 - 40, 45 - 55 keV as a template to eval- 
uate spectral variations, using the averaged response ma- 
trix weighted by individual exposures. We used a single 
power-law model in order to see possible spectral break in 
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comparison with the simple model. The values of the pho- 
ton index and flux are 2.141 ± 0.006 and 9.651^ x 10~ 9 
erg cm~ 2 s _1 (25 - 55 keV) at a 90% confidence level with 
the xVd-o.f. = 68.54/63. 

Second, we compared individual spectra with 
this model. We used the appropriate arf 

files (ae_hxd_gsohxnom_crab_20100526 . arf or 

aejixd_gsoxinom_crab_20100526 . arf ) for all 

the GSO spectra, in addition to the response 
matrices (ae_hxd_gsohxnom_20100524 . rsp or 

ae_tixd_gsoxinom_20100524 . rsp) for the aiming po- 
sitions (HXD nominal and XIS nominal). The GSO arfs 
are empirical functions to correct the normalization of 
the GSO for matching that of the PIN. implemented 
as spline functions converging to a constant above 100 
keV (Yamada et al. 2011), and the reference therein 1 . 
Figure 3 shows the ratios of individual spectra to the 
averaged best-fit power-law model. This figure clearly 
demonstrates the variations of the spectral break at ^100 
- 200 keV. For instance, the spectral slope was steep in 
the spectra of 2005.09.15, 2008.08.27, and 2012.03.14/26, 
while flat in the spectra of 2007.03.20 and 2009 - 2010. 

3.4- Quantitative Analysis of the Averaged Spectrum 

The hard X-ray spectrum (around ~100 keV) of the 
Crab Nebula is well represented by a broken power-law 
model (e.g., Yamada et al. 2011), although Jourdain & 
Roques (2009) reported that a log parabolic law (Massaro 
et al. 2000; Massaro et al. 2004) is better fitted with 
the spectrum rather than a broken power law. We thus 
tested several models for the hard X-ray spectrum using 
the HXD data. We averaged all spectra of the observa- 
tions listed in Table 1. We performed fittings with a single 
power law 

A(E) = KE r , (1) 

an exponential cutoff power law, a log parabolic law 2 , the 
Band function (Band et al. 1993), and a broken power-law 
model 

A(E) =KE r i(E<E break ), (2) 

= KE%-£ (j^j ' (E > E break ) (3) 

to the averaged spectrum. We adopted the average re- 
sponses calculated according to the exposure time of each 
observation mode (attitude and PIN epoch). Here we note 
that the spatial extension of the Crab Nebula (~ 1 arcmin 
in diameter in hard X-ray band; Pelling et al. 1987) can be 
negligible due to the HXD angular response whose field of 
view is 34' x 34' (<100 keV) and 4.5° x 4.5°(>100 keV) 
(FWHM) exhibit flat top effective areas within ~ 2 arcmin 
in diameter. 

Figure 4 shows the averaged spectrum and ratio of the 
spectrum to the best-fit models. Table 2 lists the best- 
fit parameters. The single power-law model shows large 

1 http: / /heasarc.gsfc.nasa.gov/docs/suzaku/analysis/gso_newarf.h 

2 We installed and used the logpar model as downloaded from 
http://hcasarc.gsfc.nasa.gov/xanadu/xspec/models/logpar.html 



residuals at > 100 keV to imply that a cutoff or break 
is necessary in this energy band. Although the model 
of the second smallest x 2 /d.o.f is the cutoff power-law 
model, the best-fit cutoff energy is much larger than the 
energy band of the HXD data. The broken power-law 
model, which is an empirical model for the Crab hard X- 
ray spectrum, succeeded to represent the data with the 
smallest x 2 among the five models above 100 keV range. 
We still see a residual structure in 25 - 55 keV band. This 
may suggest spectral variation in the observations. 

3.5. Quantitative Analysis of Individual Spectra 

In order to evaluate the spectral variation quantita- 
tively, we performed model fittings for individual spectra 
with the broken power-law models, which provides the 
smallest x 2 f° r the averaged spectrum (Section 3.4). As 
shown in Figure 5 and Table 3, all the spectra are repro- 
duced by the broken power law model and we also note 
that the residual structure seen in the averaged spectra 
is disappeared in these time resolved spectra. Figure 6 
shows the distribution of derived Ti and T2, which are 
photon indexes below and above the break energy, re- 
spectively. Figure 7 represents the derived break energy 
(Ebreak) from each observation. In Figure 6 and Figure 7, 
we evaluate the errors of the best fit values including sys- 
tematic errors due to the reproducibility of the GSO NXB, 
by varying the derived NXB models by ±1%. Hypotheses 
that Ti, Eb rea k and T2 were constant are safely rejected at 
>99.9%, 99.7% and 91.7% confidence levels, respectively, 
taking into account both statistic and systematic errors. 
The very conservative upper limit of Ti — T2 is 2.54 consid- 
ering the upper limit of Ti and the lower limit of r 2 at the 
90% confidence range. The upper limit of Ti — r 2 will be 
reduced < 1.29 if we exclude the data point at 2006.03.30, 
whose Ti — T2 is exceptionally large, 2.54, mainly because 
the value of T\ allows Ti = 0. We fitted Ti — T 2 with a 
constant as a function of MJD, and found that the best-fit 
constant value is 0.16 ± 0.02 at a ler error with x 2 /d.o.f. 
= 23.74/16. In Fig. 6 and Table 3, the errors increase by 
time because the level of NXB is increased due to accu- 
mulated events from radio active nuclei by cosmic-rays in 
GSO crystals. 

We also examined an exponential cutoff power-law 
model for the spectrum. Figure 8 shows the best-fit pho- 
ton indexes and cutoff energies (E cut ). We cannot obtain 
conclusive results with this model, as all best-fit values 
for E cu t fairly exceed the energy range of the HXD data 
(500 keV). 

4. Analysis of the Pulsed Component 

4-.1. Pulse Profiles 

To understand how the pulsed component of the pul- 
sar affects the variation of total emission from the Crab 
Nebula, we investigated the fluctuations of pulse profiles 
and pulsed flux of the Crab pulsar. 

2 \ We adopted the pulse period (P) and its time- 
derivative (P) from the Jodrell Bank Crab Pulsar Monthly 
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Ephemeris 3 . Epoch (phase=0) was defined as the ar- 
rival time of the first peak in the radio band. We used 
aebarycen (Terada et al. 2008) to correct photon ar- 
rival times at UTC in orbit to those at the solar sys- 
tem barycenter, assuming that the Crab pulsar is lo- 
cated in the position of R.A. = 5 h 34 m 31f97232, Decl. = 
+22°00'52"0690 (J2000), on which is the same coordinate 
used by the Jodrell Bank Crab Pulsar Monthly Ephemeris. 

Figure 9 shows the pulse profiles observed with the HXD 
at the HXD nominal pointing position shown in count s _1 . 
Here, we defined On and Off phases at phase 0.9 - 1.5 and 
0.55 - 0.85, respectively. The averaged counts at the Off 
phase are subtracted in the plot. The relative intensity 
of the first pulse peak with respect to the second one in- 
creases with the energy. The X-ray pulse peak leads to the 
first radio pulse peak (phase = 1 in Fig. 9) by ~ 300 /is, 
which is consistent with results previously reported by 
other satellites (Rots et al. 2004; Kuiper et al. 2003) 
and the initial observations with the HXD (Terada et al. 
2008). The pulse profiles were consistent with each other 
within 2.7cr of statistical errors of each phase bin through 
all observations in 25 - 300 keV. 

4-2. Pulsed Flux 

We searched for variations of the pulsed flux using 
all data at both nominal pointing positions. Count 
rates of the pulsed components, "On— Off" are plotted 
in Figure 10. Numerically, the best fit values in fitting of 
data in Figure 10 with a constant model are 1.68 ±0.01 
count s" 1 and x 2 /d.o.f. = 23.03/16 at PIN 25 - 40, 45 - 
55 keV, 6.42 ±0.04 count s" 1 and x 2 /d.o.f. = 38.96/16 
at GSO 50 - 300 keV with the errors at 90% confidence 
level. Hypotheses for the constants are rejected at 88% 
and 99% confidence levels, respectively. The pulsed com- 
ponent contributes only about 20% of the total flux, and 
thus, the possible ~ 1% of variations are hard to account 
for the total variation of the hard X-rays shown in Fig. 1 
with the fluctuation of the pulsed component. 

5. Discussion 

In Section 3.1, we demonstrated the long-term photon 
flux variation measured by the HXD, which is consistent 
with those measured by other X-ray satellites within 3 
- 4% (see also Section 3.3). The relative peak-to-peak 
amplitude of the flux variation was 6.3% in 25 - 55 keV. As 
indicated in the hardness ratio (Section 3.2), we evaluated 
the spectral variation quantitatively (Section 3.5) to find 
that the difference between the photon indexes below and 
above the break energy varies. Because the flux of the 
pulsed component contribute only ~ 20% of the total flux 
(Section 4), the variation seen in the total flux may be 
caused by the nebular component (the Crab Nebula) . The 
following sections discuss the origin of the variation of X- 
ray fluxes and energy spectra. 

3 http:/ /www.jb. man. ac.uk/~pulsar/crab. html 
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5.1. Comparison with Previous Works 

Our results obtained by the Suzaku HXD confirmed 
year-scale flux variation reported by Wilson-Hodge et al. 
(2011) using RXTE, INTEGRAL, Swift and Fermi. They 
also claimed softening in MJD 54690 - 55390 in flux de- 
creasing at <100 keV. Our results agree with their argu- 
ments regarding the general trend during this time inter- 
val, although the HXD results do not show simple soften- 
ing in detail (Fig. 2). In the range of 100 - 300 keV, 
they reported less decrease than in the softer band in 
MJD 54690 - 55390 with INTEGRAL /ISGRI. We ob- 
served similar phenomena in the HXD data, which showed 
less flux variation in the 100 - 500 keV than those of the 
25 - 55 and 50 - 100 keV bands. However, the PIN count 
rate seems to exhibit slightly larger values than those of 
other instruments after MJD 55500 by about 3 % (Fig. 1). 
This may be due to secular degradation of charge transfer 
efficiency in PIN detector. In order to examine the pos- 
sibility, we checked the PIN gains using Gd K line from 
GSO crystal and found the PIN gains in this duration 
might be reduced by —0.5 ±0.7 %. Assuming the Crab 
spectrum this corresponds to the tolerance of — 1.7j^g % 
in the photon flux, which can explain the apparent dis- 
crepancy. But we stress that a hypothesis that Ti values 
were constant is rejected at a > 99% confidence level, even 
if we omit the data after MJD 55500 during which the PIN 
gain might have been degraded. 

5.2. Origin of the Variable Component 

The X-ray radiation of the Crab Nebula is generally 
regarded as synchrotron radiation because of the non- 
thermal spectrum and the strong polarization. From our 
analysis, the averaged hard X-ray spectrum is represented 
by the broken power-law model with Ebreak ~ 100 keV 
(Section 3.4). In the PWN, electrons are re-accelerated 
and continuously injected from the termination shock, and 
then cooled down via synchrotron radiation. Assumed 
synchrotron radiation, the higher energy electrons are 
cooled in the shorter cooling time (r c ). According to the 
nomenclature of Sari et al. (1998), when r c of the mini- 
mum Lorcntz factor electrons is shorter than a duration 
(to) for each spectrum radiation, which presumably corre- 
sponds to the electron passage time through the emission 
region, the regime is called "fast cooling" regime. If it is 
the case of "fast cooling" , the observed Ebreak corresponds 
to either of v a , v c or v m , and the observed Ti is either 1, 
—2/3, or —3/2. However, none of those values arc ac- 
cepted by any of the measured Ti except for 2006.03.30 
(Tabic 3), meaning the observed spectra reject the "fast 
cooling" regime. We note that Ti of 2006.03.30 has a 
large error as Ebreak cannot be well derived. On the other 
hand, in the case of the "slow cooling" regime (to < t c ), Ti 
should be either 1, —2/3, or — (p + l)/2, and Ebreak = v a , 
v mi or j/ c , respectively. The observed values of Ti accept 
only the case of Ebreak = v c in the "slow cooling" regime. 

Although Ar = Ti — T2 should be 0.5 in the slow cool- 
ing regime since Y2 should be —(p + 2)/2, 14 out of 17 
data of ours would not accept the expected value assum- 
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ing a simple one-zone model (Ar = 0.5). However here 
we have to notice that significant spatial and temporal va- 
rations are reported in the soft X-ray band (Grcivcldinger 
& Aschenbach 1999; Hester et al. 2002; Mori et al. 2004a), 
which naturally implies multicomponent variable spectra 
also in the hard X-ray band. Moreover the filamentary 
structures seen in the Chandra images (Hester et al. 2002; 
Mori et al. 2004b) suggest that magnetic fields are differ- 
ent by locations in the Crab Nebula. On the other hand, 
the photon indexes of X-ray spectra are about 2 anywhere 
in the Crab Nebula (Mori et al. 2004b), energy spectra of 
re-accelerated electrons are thought to have almost con- 
stant indexes. Considering these facts, the multicompo- 
nent spectra are observed as an ensemble of power-law 
components and broken power-law components (with Ar 
= 0.5) in the limited observed energy band. Consequently, 
we observe variable broken power-law spectra with < Ar 
< 0.5 as superpositions of these components in the energy 
band. In fact, all of our data accept the expected range 
of < Ar < 0.5. 

Magnetic field strength B related to synchrotron radi- 
ation is given as: 



B 



100 fiG 



1 year 



-2/3 



E 



100 kcV 



-1/3 



(4) 



where E p h is the photon energy. The synchrotron bolo- 
metric intensity (Poc B 2 n(~f)~f 2 ) is assumed to be propor- 
tional to the hard X-ray flux, where 7*1(7) is the electron 
number density as a function of the Lorentz factor of elec- 
tron 7. In the slow cooling regime, r c should be compa- 
rable with the timescale of passages of electrons that are 
pertinent to v c . From our results, the timescale of the 
flux variation is expected to be a few days (Section 3.1). 
By adopting a timescale of five days, for example, we can 
estimate the value of B to be ~1.7 mG from Equation 4. 
Because the cooling timescale can be regarded as the pas- 
sage time scale, the typical scale of the variable emis- 
sion region is estimated to be around 0.35 - 1.25 light- 
days (~(3.2 - 8.1) xlO -4 pc), under the assumption of 
the speed of plasma in the PWN measured with Chandra 
(0.07 - 0.25c; Mori et al. 2004b). Although the scale of the 
variable region is ~ 0.1% of the whole X-ray emission re- 
gion, the mismatch between the occupancy (about 0.1% of 
the whole X-ray emitting region) and the factor of the flux 
variations (8% of total flux) can be solved when the vari- 
able region has ~ 10 times higher magnetic field strength 
than those in others, because P oc B 2 as mentioned before. 
Of course, because P has a dependency with 71(7) or 7, 
we may assume that a multiple component has a different 
density or Lorentz factor, but we do not discuss these pos- 
sibilities here due to limitation of information which we 
have currently. We must wait for hard X-ray images to 
conclude the entire morphology of the variable emission 
region and the magnetic fields. 
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Table 1. Suzaku Observations of the Crab Nebula 



Obs. ID 


Date 


MJD 


Exposure* 


Nominal position 


PIN epoch 


Count rate^ 




100023010 


2005/09/15 


53628 


9.9 kscc 


HXD 


1 


7.141 ± 0.028 ±0 


.003^ 


100023020 


2005/09/15 


53628 


12.4 ksec 


XIS 


1 


7.115 ± 0.026 ±0 


,004 § 


101004010 


2006/03/30 


53824 


Q 1 Vqpp 


XIS 


1 


1 . 1U^< —1— U.UOU —1— V / 


004§ 


1 01 004090 




ooo.^ y 


XZi.rJ I\.Dt;C 


XTS 


l 


J .uou m u.uzu m u 


uuo 


1 01 00S01 


ZUUU/ U^i/ uo 


ozoou 


OQ 1 L-cpr* 

Zc/.l liaLL 


1 1 . \ 1 / 


1 
1 


u.yuo in u.uiu m u 


.uuo 


1 01 01 001 


zuuo / uy / uo 


uoyoo 


ICQ Ironr 
±0.0 KSCC 


YTQ 


9 


« qqq _i_ n 091 -I- fl 


.uu^ 


1 0901 Q01 


9nn7 /n^ /9n 


o^i ( y 


An ^ l-cdp 

^U.rJ KSCC 


1 1 . \ i / 


Q 
O 


a 7qq _i_ n m 3 4- fl 
u. i oo m u.uio m u 


.uuo 


103007010 


2008/08/27 


54705 




XIS 


4 


7 040 4-0 nif! + fl 


003 


103008010 


2008/09/01 


54710 


SI Icspo 


HXD 


4 


6 966 ± 015 ± 


003 


104001010 


2009/04/02 


54923 


31.3 ksec 


HXD 


5 


6.888 ± 0.015 ±0 


.003 


104001070 

j- yj t: yj yj J. yj i \J 


2010/02/23 


55250 


1 K 1 lesec 


XIS 


8 


6 795 ± 022 ± 


.004 


105002010 


2010/04/05 


55291 


31.3 ksec 


XIS 


9 


6.794 ± 0.016 ±0 


.004 


105029010 


2011/03/21 


55641 


34.2 ksec 


XIS 


11* 


6.867± 0.015 ±0 


.004 


106012010 


2011/09/01 


55805 


33.9 ksec 


XIS 


11 


6.908 ± 0.015 ±0 


.003 


106013010 


2012/02/28 


55986 


34.0 ksec 


XIS 


lit 


6.876 ± 0.015 ±0 


.003 


106014010 


2012/03/14 


56001 


42.2 ksec 


XIS 


lit 


6.936 ± 0.013 ±0 


.003 


106015010 


2012/03/26 


56012 


22.9 ksec 


XIS 


lit 


6.961 ± 0.018 ±0 


.003 


Total 






437.7 ksec 











* Dead time corrected exposure time. 

t Count rate of the PIN in units of count s — 1 @ 25-55 keV. Errors are statistical and systematic errors of the NXB (3%) 
on la levels. Each count rate of the XIS nominal is adjusted to the HXD nominal. 

* Unusual values of the PIN epoch are selected given the special PIN LD setting during observations. 
§ Count rates of Epoch 1 and 2 are normalized to Epoch 3. 
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i 

a 

£ 



1.12 
1.1 

1, 

1.06 
1.04 
1.02 
1.00 
0. 

0.96 
0.94 
1.12 
1.1 
1.08 
1.06 
1.04 
1.02 
1.00 
0.98 
0.96 

0.94 
1.12 

1.1 

1.08 

1.06 

1.04 

1.02 

1.00 

0.98 

0.96 

0.94 




■ Sii2aku/HXD-PIN 25-55 keV 
o RXTE/PCA 1 5-50 keV 
O Swifi/BAT 14-50 keV 
□ INTEGMVIBIS 20>0 keV 
O ffermi/GBM 15-50 keV 
O INTEGM7JEMX2 1015 keV 
MAXI/GSC 10-20 keV 



■ Suzaku/HXD-GSO 50-100 keV 
O Swift/BAT50-100 keV 
□ INTEGRL/IBIS 5O00 keV 
O frnni/GBM50-100keV 




■ Suzaku/HXD-GSO 100-500 keV 
::: INTEGRUBIS 100-300 keV 
O Fcrmi/GBM 100-300 keV 



- 1 





51000 52000 53000 54000 

MJD 



55000 



56000 



Fig. 1. Long-term light curves with the Suzaku/HXD in 25-55 keV(top), 50-100 keV (middle), and 100-500 keV (bottom) overlaid 
on data of Swift, RXTE, INTEGRAL, Fermi and MAXI. The data obtained by RXTE, Swift, INTEGRAL and Fermi before MJD 
55000 arc from Wilson-Hodge et al. (2011). Swift after MJD 55000 and MAXI (Matsuoka et al. 2009) are taken from their archival 
data [(Swift) http://swift.gsfc.nasa.gov/docs/swift/results/transients/, (MAXI) http://maxi.rikcn.jp/]. The error bars of the HXD 
include lcr statistical errors and la systematic errors of NXB (Fukazawa et al. 2009). Each flux is normalized to mean flux in the 
time interval of MJD 55000-55500. 
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6.4 6.5 6.6 6.7 6.8 
count/sec (25-40, 45-55 ke V) 



0.68 
0.67 
0.66 
0.65 




20.2 204 20.6 20.8 21 21.2 21.4 21.6 21.S 
count/sec (50-500 kcV) 



Fig. 2. Count rates vs. hardness ratios of PIN 30-40, 45-50 keV/25-30 keV (left) and GSO 100-500 keV/50-100 kcV (right). The 
horizontal axes show count rates and the vertical axes show hardness ratios of count rates. Data marked with open circles in the 
left graph are only references because different bias voltages of the PIN were set in these observations. The error bars represent 
statistical errors (lc). 



Table 2. Best-fit models of the averaged spectrum 



Model name 


Parameter name 


Value (error) 


Single power law 


r 

Flux/ 
Flux/ 
X 2 / d.o.f. 


-2.125 ±0.001 
0.970 ±0.009 
7 865 +a008 

1 -000-0.007 

797.16 / 87 


Cutoff power law 


r 

E cut (keV) 
Flux/ 
Flux/ 
X 2 / d.o.f. 


-2.075 ±0.004 
1396±^| 

n Q79S+0-0008 
u-»'^o_o .0010 
7 OKI +0-008 

431.04 / 86 


Log parabolic law 


a 

P 

pivot E (kcV) 
Flux/ 
Flux/ 
X 2 / d.o.f. 


-2.150 ±0.003 
-0.059 ±0.006 
20 (fix) 
q73q+ 0010 

u.»io»_ 0009 

7.948 ±0.011 
540.28 / 86 


Band function* 


a 

P 

E peak (keV) 
X 2 / d.o.f. 


Z.U(O_ 004 
> -2.387 
1378jjf 7 
429.25 / 85 



Broken power law 



Ebreak (keV) 

r 2 

Flux/ 
Flux/ 
X 2 / d.o.f. 



-2.115±0.002 

1341L 
-2.26 ±0.02 

q700 +0 0009 
u.a/uu_ .ooo8 

7.916 ±0.009 
311.14 / 85 



All errors arc presented at 90% confidence levels. 
* in units of xl0 — 8 erg cm~ 2 s" 1 at 25-55 keV 
t in units of xlO -9 erg cm -2 s _1 at 50-100 keV 

t Flux of the Band function cannot be derived as /3 and E p< , ak have too large 
errors. 
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20050915 (HXD nom) 



-V*-t- 







-+- 



20060404 (XIS nom) 



' + +++f* 



-+- 



20060405 (HXD nom) 



^ 



20060905 (XIS nom) 



20070320 (HXD nom) 



20080827 (XIS nom) 



-H 



20080901 (HXD nom) 
I. 



. 20090401 (HXD nom) 
.1 



1 20100223 (XIS nom) 



. 20100405 (XIS nom) 



20110321 (XIS nom) 
Jl 



WW 



20110901 (XIS nom) 



1 201 20228 (XIS nom) 



20120314 (XIS nom) 



■_ 20120326 (XIS nom) 



300 50 

Energy (keV) 
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Fig. 3. Ratios of background-subtracted spectra of individual observations to the best-fit single power-law model of the averaged 
PIN spectrum (25-40, 45-55 keV) of all observations. The error bars represent statistical errors (1<t). 
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Fig. 4. The top panel shows the averaged spectrum through 2005-2012 with the Suzaku HXD. The solid line represents the best-fit 
single power-law model and crosses are data with lcr statistical errors. Below five panels show the ratios of the data to best-fit 
models of a single power law (power), an exponential-cutoff power law (cutoffpl), a log parabolic law (logpar), Band function (grbm), 
and a broken power law (bknpower). 
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Fig. 5. Ratios of data to best-fit broken power-law models of individual spectra 
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Table 3. Best-fit broken power-law model of individual spectra 



[Vol. , 



Flux/J 

hi). 06 



2~T 



ObsID 



Date 



MJD Ti 



_9 nQ o+0. 014+0 
z.u»o_ oio-o 

_9 nQ^+ 020 +° 

z.u»O_ .015-0 
-9 1OC.+0. 018+2 
z - lzd -0. 010-0 

-2 inq+0-008+o 

Z.1U»_ 008-0 

-9 nQ4+°- 005 +° 
z.u»1_ .005-0 

-2 iofi+ 012 +° 

Z.1UU_ oio-O 

-9 nQ*?+°- 007 + 

z.u»O_ 007-0 
_9 1 9H+0-006+0 
Z.1ZU_ . 006-0 

_9 in^+°- 005 + 

Z.1UO_ 005-0 

_9 im+°- 007 + 
z - luJ -0. 007-0 
_ 9 -i 90+O.OII+O 
Z.1ZO_ 008-O 
-9 1 9/1+0-006+0 
z.izi_ 006-O 
_9 1 QO+0. 005+0 
z - loo -0. 005-0 
_9 1 1 7+0-007+0 
z - 11 ' -0.007-0 
— 9 1 70+0. 023+0 
'°-0. 035-0 
-9 1 9o+°- 005+0 
z ' 1 ^ d -0. 005-0 
_9 1 Qo+0. 007+0 
Z - LOO -0.007-0 



TTGT 

004 
001 
005 
107 
005 
000 
005 
011 
004 
000 
004 
003 
007 
01 IS 
00!) 
013 
007 
000 
006 
000 
064 
000 
159 
011 
005 
005 
006 
034 
000 
DOS 
000 
009 
010 



E break (keV) 
1 ori + 2i! + 5 

lzu -25-10 
09+29+8 
yz -20-6 
1 09+38+29 
ioz -56-76 



r 2 

RT 
'-0 
HO 
-0 
^+0 

-2.351° 
-2.28+° 
-2.18 



Flux, 



.08+0.04 

10- 0.10 
.03+0.03 
06-0.08 
09+0.03 

14- 0.14 
.09+0.00 

11- 0.12 
05+0.09 
06-0.09 
.03+4.96 
03-0.05 
02+0.04 

03- 0.10 
.13+0.15 
21-0.40 
12+0.06 

08- 0.14 
.03+0.00 

04- 0.07 
13+0.06 
79-0.38 

.06+0.12 

09- 0.19 
35+0.20 
.22-0.57 
.03+0.06 

05- 0.07 
01 + 1.15 
01-0.28 
.11+0.13 

15- 0.28 
18+0.16 
25-0.34 



1 oo5 +tU)Ub 

l.UUO_o.oo6 

1 000+ 006 
i.uuu_ 006 

1 oorr - 006 

i - uuu +0.006 

n nqq+0.005 
u.a»o_ 005 
n QS1 +0.003 

U.»O±_ 003 

n Q82+ 004 

U.»OZ_ 005 

0.949iH°o! 

0.987j« 
n 070+0.003 

u - y ' y -0.003 
ORO+0.003 
u - aua -0.003 

o.948iS:88» 

0.9501°;°$ 

0.9581S 

0.965l°;°° 3 3 
954+0.0O6 

070+0 003 

\J.v u — 0.003 

n 070+0.004 

u - b " u -0.004 



X 



100023010 
100023020 
101004010 
101004020 
101003010 
101010010 
102019010 
103007010 
103008010 
104001010 
104001070 
105002010 
105029010 
106012010 
106013010 
106014010 
106015010 



2005/09/15 
2005/09/15 
2006/03/30 
2006/04/04 
2006/04/05 
2006/09/05 
2007/03/20 
2008/08/27 
2008/09/01 
2009/04/02 
2010/02/23 
2010/04/05 
2011/03/21 
2011/09/01 
2012/02/28 
2012/03/14 
2012/03/26 



52628 
53628 
53824 
53829 
52830 
53983 
54179 
54705 
54710 
54922 
55250 
55291 
55641 
55805 
55986 
56001 
56012 



14F 



Values of the first errors of Ti , Et rea k and 
the sum of statistical and systematic errors 



I +16+1 
-19-0 
1 oc + 12+1 
i00 -15-34 
Q9+17+304 
az -17-0 
qq+20+12 
Jy -ll-7 

166±il±l 

1 1 1 +13+0 
111 -ll-0 
1 yq+171+0 
1 ' ,3 -64-76 
145+27+0 
1 ^ :J -26-89 
947+30+0 
-19-143 

iU8 -15-93 
^r + 10+242 
^ J -12-0 
i 66 +17+6 
iuu -15-6 

173±i»+ 4 



-2.30] 
-2.21; 
-2.26^ 




+0 

-2.19±§ 

-2.601°, 

-2.32+° 

-2.251° 

-2.261°, 
^+o 

-1.44±§ 

-2.201° 
1+0 

-2.54+° 
1+0 

-0 



-2.26 H 



-2.11" 



-2.74 H 



8.28^ 



.03^ 



8-32_o.o6 

2+0.05 
-0.1 
o 1 1 +0.06 
8-ll_0.07 

o -1 O+0.05 
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Fig. 6. Ti vs. T2 plots of the best-fit broken power-law models for individual spectra. The small error bars (bold lines) represent 
statistical errors and the large error bars (fine lines) represent systematic errors of the NXB uncertainty (see text). All errors are on 
a 90% confidence level. 
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Fig. 7. Time (MJD) vs. -Ej, rea fe of the best-fit broken power-law models for individual spectra. The error bars represent statistical 
errors plus systematic errors of the NXB uncertainty (see text). All errors are on a la level. The solid line represents the best-fit 
constant function. 
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Fig. 8. r vs. E cu t of the best-fit cutoff power-law models for individual spectra. The dashed line shows E cu t = 500 keV. The small 
error bars (bold lines) represent statistical errors; the large error bars (fine lines) represent systematic errors of the NXB uncertainty 
(see text). All errors are on a 90% confidence level. 
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Fig. 9. Pulse profiles of the Crab pulsar observed with the HXD nominal position. Vertical axes show the count rates with 
subtraction of the average Off phase values. The error bars represent statistical errors (lcr). 
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Fig. 10. Count rates of the Pulsed components of the Crab pulsar in 25—40, 45—55 keV of the PIN (top) and in 50-300 keV of the 
GSO (bottom). The dashed line in each panel represents the best- fit constant function in each energy band. All errors arc on a 90% 
confidence level. 



